固有多項式 $p(\lambda) = \det(A - \lambda I)$ は固有値を定義する理論的基盤ではあるが、高次元システムでは数値的に「不安定」で計算上も非効率である。実用的な応用（例えば波動伝播のストゥルム・リウビル問題）において、多項式の根が係数の摂動に極めて敏感なため、直接展開するよりも他の選択肢が優先される。
連続波から離散行列へ
弦や膜の振動は次の波動方程式によって支配される：
$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$
解 $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$ を求めるには、 ストゥルム・リウビル系を解く必要がある。
$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$
離散化の複雑さ
演算子の離散化により、$Aw = -0.04 \frac{\rho}{p} \lambda w$ のような行列方程式が得られる。$4 \times 4$ の三重対角行列では $p(\lambda)$ は扱いやすいが、メッシュが細かくなる（$n$ が増加する）と二つの壁に直面する：
- アーベル・ルフィニの制限： 次数 $n \ge 5$ の多項式の根に対して代数的解は存在しない。
- 丸めの感度： 高次元システムでは、行列の要素の $10^{-10}$ 桁のわずかな変化が、固有値を桁違いに変化させる（ウィルキンソンの多項式現象）。
数値的必要性と専門ライブラリ
専門的な数値ライブラリ（IMSL、NAG）は原始的な固有多項式を避ける。代わりに、近似のために反復的手法を使用する：
- IMSL ライブラリ： 線形最小二乗法、立方スプライン、高速フーリエ変換を使用する。
- NAG ライブラリ： 最小二乗多項式近似および $l_1/l_{\infty}$ 意味でのフィットを採用する。
系統 $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$ における固有値の近似では、根の探査ではなく離散最小二乗法と反復的発見に頼る。
🎯 理論的ツールと数値的危険性
固有多項式 $p(\lambda) = \det(A - \lambda I)$ は証明において重要だが、計算では危険である。物理における実用的な固有値問題は、安定性を保つ反復変換（例：QR法）によって解決される。